Transition metal ad-atoms on graphene: Influence of local Coulomb interactions on 

chemical bonding and magnetic moments 
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We address the interaction of graphene with 3d transition metal adatonis and the formation 
of localized magnetic moments by means of first-principles calculations. By comparing calculations 
within the generalized gradient approximation (GGA) to GGA+U we find that the electronic config- 
uration and the adsorption geometries can be very sensitive to effects of local Coulomb interactions 
U in the transition metal d-orbitals. We find high-spin configurations being favorable for Cr and 
Mn adatoms independent of the functional. For Fe, Co and Ni different electronic configurations are 
realized depending on the value of the local Coulomb interaction strength U. Chemical control over 
the spin of the adatoms by hydrogenation is demonstrated: NiH and CoH adsorbed to graphene 
exhibit spin 5* — 1/2 and 5=1, respectively. 



I. INTRODUCTION 

Graphene — a single layer of C-atoms arranged to a 
honeycomb lattice — has developed from a theoretical 
model systeroii^ to experimental reality^. It is stabilized 
by a network of sp^ cr-bonds and conjugated tt bonds, 
which make this material chemically inert. Its low energy 
electronic structure is determined by two tt bands inter- 
secting linearly in the corners of the Brillouin zone. The 
Fermi level of the pristine material lies at these cross- 
ing points — called Dirac points. As a consequence, 
graphene exhibits Dirac fermion excitations and a lin- 
early vanishing density of states (DOS) at low energies. 
Carrier concentrations have been proven tunable — both 
for electrons and holes — by more than an order of mag- 
nitude by chemical^"— as well as electrostatic doping^. 
It can be anticipated that this control over the electronic 
density (of states) implies a general tunability of electron 
correlation phenomena in graphene. 

Adsorbates carrying a magnetic moment, like 3d tran- 
sition metal adatoms, and binding to graphene present 
a chemical route to localized magnetic moments in 
graphene and possibly allow us probing tunable many 
body physics in this material. Various studies have ad- 
dressed the adsorption of 3d adatoms by means of density 
functional (DFT) calculations^"— . 

Many Sd-adatoms including Fe, Co and Ni were pre- 
dicted to adsorb in the center of a graphene hexagon 
(h-site)^i^. At h-sites, the coordination number is max- 
imized for the adatom. Its d-orbitals with non-zero z- 
component of the orbital angular momentum {l^ ^ 0) 
couple to the low energy graphene states, whereas the 
s-orbital of the adatoms remains decoupled from the 
low energy states a this sitei^ii^. Away from the h- 
sitc position, the s-orbital can hybridize with low energy 
graphene states but forming a cr-bond with one of the car- 
bon atoms requires energy for breaking graphene's sp^- 
network2£. Moreover, the reduced coordination decreases 
the d-electron contribution to bonding. The adsorption 
site of adatoms relies hence on balancing the contribu- 
tions to chemical bonding from the s- against those from 



the d-orbitals. The spatial extent of the latter orbitals 
is strongly influenced by the local Coulomb interaction. 
Therefore, the energy gain upon bonding can be very 
sensitive to effects of on-site Coulomb repulsion. 

In this paper, we investigate how the interaction of 3d 
transition metal adatoms with graphene is affected by lo- 
cal Coulomb interactions and consider Cr, Mn, Fe, Co, 
and Ni atoms adsorbed to graphene. In section |ll] de- 
tails of our DFT calculations are given. In section IIIIl 
we present our results on adsorption geometries and the 
electronic structure of 3d adatoms on graphene. The in- 
fluence of local Coulomb interactions is discussed. In 
section IIV[ we show that the magnetic state of adatoms 
on graphene is very sensitive to their chemical environ- 
ment by considering hydrogenated Co and Ni adatoms 
on graphene. Conclusions are given in section IVl 



II. COMPUTATIONAL METHOD 

For an ab-initio description of the transition metal 
adatoms on graphene we performed DFT calculations on 
3x3 and 4x4 graphene supercells containing one adatom 
using the Vienna Ab Initio Simulation Package (VASP)^ 
with the projector augmented wave (PAW) ^^'^'^ basis 
sets. To judge the role of local Coulomb interactions 
in the 3d orbitals of the adatoms, we employed a gen- 
eralized gradient approximation functional (GGA)^ as 
well as GGA+U with [/ = 4eV, J = 0.9 eV, which is 
typical value for 3d-transition metal impurities in metal- 
lic host o^^i^^ . In addition, LDA calculations were per- 
formed. The calculations were done with plane wave cut- 
offs 800 eV and 928 eV for the GGA and LDA functionals, 
respectively. We obtained fully relaxed structures for all 
of these functionals with the forces acting on each atom 
being less than 0.02 eVA 

The binding energies Eb reported in the following were 
obtained using the standard formula Eh = Eg + i?ads — 
-E-g-ads, where Sg-ads, Eg and i?ads are, respectively, the 
total energies of the graphene-adsorbate system, a pris- 
tine graphene sheet and the isolated adsorbate. As iso- 



latcd atoms with open electronic shells arc strongly cor- 
related systems, the GGA or GGA+U values for i?ads can 
be unreliable and it is hence most meaningful to consider 
binding energy differences where i?ads cancels out. 



III. MAGNETISM AND STRUCTURE OF 3rf 
ADATOMS ON GRAPHENE 

We recapitulate the example of Co which has been sub- 
ject of prior studies^ii^iiiii^iii and then generalize to Cr, 
Mn, Fc and Ni adatoms on graphene. 

To judge which adsorption geometries are possible, we 
firstly relaxed Co adatoms adsorbed to high-symmetry 
sites the graphene lattice: Co above the center of a 
hexagon (h-site), Co on top of a C-atom (t-site) and 
above the middle of the bridge (b-site). Using a 4 x 
4 graphene supercell, 100 meV first order Methfessel- 
Paxton smearing and a 6 x 6 k-mesh for sampling the su- 
percell Brillouin zone, we calculated to total-energy and 
the total magnetic moment of the supercell for the Co 
atom being moved along the path connecting the relaxed 
adsorption geometries^l. As can be seen from Fig. [1] 
LDA as well as GGA predict Co adsorbing to a h-site 
with t-site or b-site adsorption being > 0.5 eV higher in 
energy. In all cases, the supercell magnetic moment is 
near tti « 1 /i^ suggesting an electronic configuration of 
Co close to spin S — 1/2. These results arc in agreement 
with prior studies^ii^ii^. 

Orbital dependent on-site Coulomb repulsion, as ac- 
counted for in GGA-I-U in a mean field manner, affects 
the predicted configuration of Co on graphene drastically: 
With U = 4eV, we find that the potential energy land- 
scape is reversed as compared to the GGA / LDA pre- 
diction (see Fig. [T] middle). For U = 4eV, J = 0.9eV 
the global minimum energy is found for Co in a high spin 
state with S ~ 3/2 on a t-site, which is 0.2 cV and 0.08 eV 
lower in energy than the h- and b-sites, respectively. 

The local density of states depicted in Fig. [5] show 
that the electronic configuration of Co markedly changes 
upon increasing the local Coulomb interaction U. In the 
GGA low spin solution the Co 4s-orbital is unoccupied 
and approximately one hole resides in the Co d-orbitals 
with El symmetry {dxz, dy^), i.e. Co is close to a 3d^4s° 
configuration. In contrast, at C/ = 4eV Co appears in a 
3d^4s^ configuration. 

This change in the electronic properties comes along 
with a change in equilibrium heights of the Co adatoms: 
With the increase of the Co spin from ^ 1/2 to ^ 3/2, the 
equilibrium height of Co above the sheet increases from 
1.5 A to 1.8 A for the adatoms at an h-site as well as 
from 1.8 A to 2.2 A for the adatoms placed at the t-site. 
This observation is in agreement with Ref. Il6l . where 
the B3LYP hybrid functional has been employed. The 
change in geometry can be traced back to the more ex- 
tended Co 4s-orbital becoming occupied, when increas- 
ing the Hubbard-U and the Co spin. Balancing electrons 
between the Co 3d and 4s-orbitals is hence decisive for 



determining the electronic properties and the adsorption 
geometry of Co on graphene. 

This balancing of electrons between the 3d and 4s- 
orbitals is related with the intcrconfiguration energies 
like AEic = £;(4si3<i"~i) - EiAs^Sd"--^). For free 3d- 
atoms it is known that semilocal functionals like LDA or 
GGA tend to underestimate AE'ic and generally favor 4s 
occupancies lower than experimentally observed^^. For 
Co adatoms on graphene, it is not yet clear whether the 
LDA/GGA configuration with Co close to a 3d^As'^ or 
the GGA+U 3(i^4s^ configuration is closer to the exper- 
iment or whether even both configurations might be re- 
alized. The example of Co adatoms, however, demon- 
strates that for a general understanding of 3d transi- 
tion metal adatoms on graphene comparing semilocal 
LDA/GGA type functionals to GGA-I-U can be useful. 

To this end, we performed GGA as well as GGA+U 
(t/ = 4eV, J = 0.9 eV) calculations for Cr, Mn, Fe, Co, 
and Ni on graphene (one adatom per 3x3 graphene su- 
percell), which allowed us to obtain fully relaxed geome- 
tries, total-energies, and magnetic properties. The re- 
sults are summarized in Table HI Two classes of adsorbed 
electronic configurations are obtained: Those with the 
supercell magnetic moment close to the free atom mo- 
ment (high-spin configuration, colored black) and those 
exhibiting a by ~ 2hb smaller magnetic moment than the 
free atoms (low-spin configuration colored red (grey)). 
The latter solutions become increasingly favorable for the 
later 3d elements: In agreement with Refs. [^, M, [llHlJ 
our GGA calculations yield Fe in low-spin configuration 
at an h-site, as well as Co and Ni in low-spin configura- 
tion at all sites. 

In all cases with stable low-spin configuration, h-site 
adsorption is favored. The LDOS (Fig. [2|) shows that in 
all low-spin configurations 4s-orbitals of the adatoms are 
unoccupied with electrons being transferred to the 3d- 
orbitals. In the low spin h-site configurations obtained 
for Fe, Co, and Ni the dxz- and dj^^-orbitals {Ei type 
orbitals) are the 3(i-orbitals with the highest energy. For 
Fe (Co), GGA predicts a spin of 5 = 1 (S* = 1/2) due 
to two (one) holes in the Ei type orbitals, which are 
potentially most important for spin fluctuations in both 
cases. 

As the example of isolated Co adatoms demonstrated, 
local orbital-dependent Coulomb interaction can strongly 
affect the balance of 3d and 4s electrons and consequently 
change the adsorption geometry as well as the magnetic 
state of the adatom. As in the case of Co, the low-spin 
solution of Fe at an h-site becomes unstable for [/ = 4 eV. 
At U = 4eV there exist only high-spin solutions for Cr, 
Mn, Fe, and Co. 

Only Ni remains in the low-spin configuration also 
for U = 4eV. There are however metastable GGA+U 
high-spin solutions also for Ni with magnetic moments 
m = 1.9 fiB and m = 1.6 fiB for Ni at b- and h-sites, 
respectively. These high-spin solutions are, respectively, 
0.27eV and 0.35 cV higher in energy than the low-spin 
h-site configuration. At U = 6eV, the b-site high spin 
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FIG. 1. (Color online) Left: Schematic illustration high symmetry adsorption sites for Co adatoms on graphene (red) and paths 
connecting these sites (green). Middle: Total energy, E, of graphene with a Co adatom as function of the Co position along 
the path depicted in the left panel. Energies obtained with GGA, LDA, and GGA-I-U (\J — 4eV, J = 0.9 eV) are shown. The 
energy for Co at an h-site is defined as E — 0. Right: Total magnetic moment of the graphene supercell with adsorbed Co as 
function of the Co position. 



Magnetic moment (hb) 




GGA 


GGA+U 




t 


b 


h 


t 


b 


h 


Cr 
Mn 
Fe 
Co 

Ni 


5.9 
5.1 
4.1 
1.1 
0.0 


5.8 
5.1 
4.1 
1.1 
0.0 


5.6 
5.2 
2.0 
1.1 
0.0 


5.9 
5.0 
3.9 
2.9 
0.0 


5.9 
5.0 
3.9 
2.9 
0.0 


5.9 
5.0 
3.6 
2.6 
0.0 



Height above sheet (A) 




GGA 


GGA+U 




t 


b 


h 


t 


b 


h 


Cr 

Mn 

Fe 

Co 

Ni 


2.4 
2.2 
2.2 
1.9 
1.9 


2.3 
2.2 
2.2 
1.8 
1.8 


2.1 
2.1 
1.5 
1.5 
1.6 


2.6 
3.7 
2.2 
2.2 
1.8 


2.6 
3.8 
2.3 
2.2 
1.8 


2.9 
3.8 
2.0 
1.8 
1.5 



Binding energy (eV) 




GGA 


GGA+U 




t 


b 


h 


t 


b 


h 


Cr 
Mn 
Fe 
Co 

Ni 


0.49 
0.35 
0.24 
0.75 
1.44 


0.50 
0.34 
0.22 
0.88 
1.52 


0.47 
0.42 
0.71 
1.35 
1.72 


0.34 
0.30 
0.29 
0.62 

0.74 


0.34 
0.28 
0.25 
0.56 
0.82 


0.33 
0.28 
0.27 
0.41 
0.97 



TABLE L Magnetic moment of supercell, height above graphene sheet and energy relative h-site for 3d transition metal adatoms 
as obtained from GGA and GGA+U with U = 4eV and J — 0.9 eV. High-spin solutions are colored back, low-spin solution 
red (grey). 



configuration becomes energetically most favorable. 

The LDOS at the adatoms and the total DOS de- 
picted in Fig. [2] illustrate bonding mechanisms depend- 
ing on the electronic configurations of the adsorbates. 
For all species in low-spin configuration there is strong 
hybridization between the d-orbitals and the graphene 
bands leading to the supercell Dirac point to being no 
more recognizable as discussed for Fe in Ref. [^ 

For Cr, both, GGA and GGA+U predict high-spin 
states with similar LDOS near the Fermi level {E = 0) 
and the Dirac point being visible for the minority spin 
spin states at w — 0.65eV and — 0.52cV, respectively. 
Given the 3x3 supercell size this indicates ionic bond- 
ing with Cr donating on the order of 0.1 — 0.2 electrons 
to the graphene bandsi^ Similarly, GGA predicts ionic 
bonding of Mn on graphene with the Dirac point shift 
indicating 0.3 — 0.5 electrons being transferred to the 
graphene bands. GGA+U yields Mn close to the atomic 
3dP4s^ electronic configuration and much weaker bonding 
to graphene. For Fe, various resonances in the vicinity of 
the Fermi level hinder to distinguish between graphene 
bands and adatom induced states making an estimate 
of doping due to Fe ambiguous. The supercell magnetic 
moment m = 3.9 ^b in the high-spin state suggests Fe 
acting as donor with « O.le being transfcred to graphene. 

The LDOS of the 4s-orbitals of the adatoms consist 
mainly of a narrow peak and indicate at 4s-electrons 



remaining localized upon adsorption of the atoms on 
graphene. In many high-spin solutions (Cr (GGA and 
GGA+U), Mn (GGA), Fe (GGA+U), Co (GGA+U), and 
Ni (GGA+U)), the 4s resonance of the adatoms is oc- 
cupied by (almost) one electron and close to the Fermi 
level. In these cases, the impurity 4s-orbitals contribute 
a magnetic moment of Ifis- Hence, when discussing the 
low energy physics (e.g. the Kondo effect) of adatoms in 
high-spin configurations on graphene, the 4s-orbital has 
to be taken into account, as it is potentially important 
for charge or spin-fiuctuations of the magnetic adatom. 

The coexistence of high- and low-spin solutions for 
Fe, Co, and Ni suggests that the spin of transition 
metal adatoms can be tunable. In the following sec- 
tion, we demonstrate that hydrogenation presents chem- 
ical means to control the magnetic properties of adatoms 
on graphene with the examples of Co and Ni. 



IV. HYDROGENATED CO AND NI ADATOMS 

Transition metal s-orbitals are, in general, highly sus- 
ceptible to chemical bonding. In bulk 3(i-elements, 
the 4s-orbital is completely delocalized causing metallic 
bonding and gives rise to a conduction electron band- 
width on the order of some 10 eV. Even in the case of 
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FIG. 2. (Color online) Local orbitally and spin resolved density of states at the adatoms (Cr, Mn, Fe, Co, and Ni from top 
to bottom) and total density of states of the supercells obtain with GGA (left column) as well as GGA+U (right column, 
U = 4eV, J = 0.9 eV). Spin-up is plotted on the positive ordinate, spin-down on the negative ordinate. The GGA-I-U high 
spin solution show for Ni is metastable at i!7 = 4 eV. All other LDOS shown here belong to structures with lowest total energy. 




FIG. 3. (Color online) Relaxed adsorption geometries of CoH 
on graphene as predicted by GGA (left) and GGA+U (right). 
For NiH, both, GGA and GGA+U yield h-site adsorption 
with virtually the same structure as obtained for CoH in GGA 

(left). 



reduced coordination, e.g. for Co adatoms on metal sur- 
faces like Cu (111) (see e.g. [30), the characteristic energy 
band width of the 4s-orbital is on the order of 10 eV. Sim- 
ilarly, for transition metal adatoms on insulating surfaces 
like CuN)^ covalcnt bonds cause distribution of 4s spec- 
tral weight over an energy range on the order of 10 eV. 

For adatoms on graphene, the situation turned out to 
be very different (see section UlI)) . where the 4s-orbital 
forms no strong covalent bond with a carbon atom un- 
derneath. This "frustration" of the transition metal 
4s-orbital can be resolved by chemical bonds to other 
species. Hydrogen is available in many experimental sit- 
uations and consequently particularly important. Here, 
we consider the interaction of H with Co and Ni adatoms, 
where the 4s-orbital is predicted by GGA to be unoccu- 
pied. 

Our structural relaxations show that H forms a chem- 
ical bond with the Ni and Co adatoms leading to a 
CoH bond length of 1.54 A (1.55 A) and a NiH bond 
length of 1.49 A (1.51 A) in GGA (GGA+U; U acting 
on transition metal d-orbitals; U = 4eV, J = 0.9 eV). 
The fully relaxed adsorption geometries are shown in 
Fig. El For NiH both, GGA and GGA+U yield h-site 
adsorption, where the Cq^ symmetry is broken by the 
NiH bond. For hydrogenated Co, we obtain h-sitc ad- 
sorption in GGA but GGA+U predicting t-site adsorp- 
tion. For CoH, the adsorption heights of the transition 
metal atoms of 1.7 A (at h-sitc, GGA) and 2.1 A (t-site, 
GGA+U) are very close to the heights of Co in the re- 
spective high spin GGA+U structures. This gives a first 
indication that extended Co 4s derived states become 
partially occupied upon hydrogenation. Similarly, for 
NiH the metal atom heights of 1.7 A (1.8 A) above the h- 
sitc obtained from GGA (GGA+U) increase compared to 
the heights of isolated Ni atoms in low-spin configuration 
on graphene. Independently of the employed functional 
(GGA or GGA+U), the supercell magnetic moments are 
2.0 /iB and 1.0 /^s for CoH and NiH on graphene, respec- 
tively. 

Well in line with the structural and magnetic changes 
upon hydrogenation, the LDOS of the Co Zd electrons in 
CoH is very similar to the high-spin solutions obtained 
for Co on graphene in GGA+U. The h-site configuration 
obtained for CoH in GGA has a d-electron LDOS which 
is very similar to the high-spin h-site solution for isolated 
Co on graphene (c.f. Ref. [Ij). The t-site GGA+U LDOS 



of CoH is similar to the LDOS of isolated Co at a t-site 
(Fig. [5|). The major difference between hydrogenated 
and isolated Co on graphene is that there are no mainly 
4s-derived states pinned to the Fermi level in the former 
case. 

For NiH, the LDOS (Fig. ^ shows that its mag- 
netic moment is due to a hole in an orbital with large 
Ni dxz , dyz-character which would correspond to orbitals 
transforming according to Ei in a fully Cg^ symmetric 
environment. Here, C^v symmetry is broken by the NiH 
bond (Fig. [3|) . Analyzing the occupation matrix of the 
Ni 3d-subspace as defined by the PAW basis sets we find 
9 eigenvalues > 0.84, i.e. 9 nearly filled Ni d-orbitals. 
One significantly lower eigenvalue of 0.42 (GGA) as well 
as 0.22 (GGA+U) corroborates that the hole leading to 
the supercell magnetic moment has largely Ni Sd char- 
acter. The projection Vei of the eigenvector j^I'o) be- 
longing to this eigenvalue onto the Ni dxzidyz-snhsp&ce 
yields {"H q\V ei\^ o) = 0.96 in GGA and 0.97 in GGA+U. 
The hole in the Ni 2>d derived states has consequently the 
major contribution from Ei type orbitals. Hence, if NiH 
on graphene gives rise to a Kondo effect, the the orbitally 
controlled Kondo effect discussed for Co on graphene in 
Ref. [ij likely generalizes the case at hand. 



V. CONCLUSIONS 

Based on first-principles calculations we discussed the 
interaction of graphene with transition metal adatoms. 
On-site Coulomb interactions can decisively affect in- 
terorbital charge transfers and thus control the electronic 
configuration, the adsorption geometry and the stability, 
particularly migration barriers, of adatoms on graphene. 

Ad-atoms with high-spin ground state and (partially) 
occupied 4s-orbital exhibit, in general, migration bar- 
riers on the order of < 0.06 eV. In contrast, adatoms 
with low-spin ground states and mainly empty 4s-orbital 
are stabilized by covalent bonds of the Sd-electrons to 
graphene and have significantly larger migration barriers 
of 0.15 — 0.5 eV. The thermal stability of adatoms can 
be probed, e.g., by scanning tunneling microscopy at dif- 
ferent temperatures and will allow to draw conclusions 
on the electronic configuration of the adatoms. This is 
particularly interesting for the cases of Fe, Co, and Ni 
where GGA and GGA+U predict different ground state 
configurations. 

The coexistence of high-spin and low spin solutions 
for Fe, Co, and Ni can be traced back to the fact that 
the 4s-orbital of the adatom is usually chemically active 
but "frustrated" on graphene: For Sd-transition metal 
adatoms on the normal metal surfaces the 4s-electrons 
arc fully dclocalized and basically part of the conduction 
electron sea, whearas for adatoms on graphene a half- 
filled 4s-orbital will be part of the impurity spin. 

Hydrogentation resolves the frustration of the 4s- 
clectrons and allows to tune magnetism of transition 
metal atoms on graphene chemically: NiH exhibits spin 
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FIG. 4. (Color online) Local orbitally and spin resolved density of states for the hydrogenated adatoms CoH and NiH on 
graphene (from top to bottom) obtained with GGA (left column) as well as GGA+U (right column, U — 4eV, J — 0.9 eV). 
Spin-up is plotted on the positive ordinate, spin-down on the negative ordinate. 



S = 1/2 on graphene, while isolated Ni adatoms are non- 
magnetic in GGA. Whether high- or low spin states are 
realized for a given transition metal adatom on graphene 
can be expected to depend very generally on the chemical 
environment. This issue including the influence of differ- 
ent substrates, e.g. insulators like Si02^i^ or metals like 
Ir— i2i, deserves future investigations. 
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